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Ai^SiRACT 

I ';-.' .-'.mal'-si-alv' horizontal U i.st btructure ci' detailed v/lr.d 
o.-oliie ,-i along Uic vortical in the liobt 20 km of tlic atmosplie re is a 
■. o rt •oally nonhomc},enoout; process. A linear stocnastic moce. is 
c'.ovoloooci lor thiC process based on tnc process covariance iuncticn. 
Ihii model is formulated through the use of a scaling hypothesis 
wr.Lch transforms the nonhomogeneous gust process into a nondimen- 
si mal gust process which is homogoneou'. in a nondimensional height 
^oordinme. The velocity scaling parameter for the gust process is 
tin: gust standard deviation, and the length scale used to nondimen- 
sionalize the altitude is the vortical space lag associated with the 
.hrst aero of the gust covariance function. State space theory is used 
to derive a digital filter from the model, which can be readily used to 
simulate gusts for space vehicle design applications. 





i .s I j cC V ! j 1.^ 1 1 •> ' .■ 

I'.uir'i ht r ik ii'i:' v iiaiMv u-rl::i ic »>r ;• vertical profile- of the 
.o rir.onial wind in the first 2<) km of the atmosphere is its stochastic 
■. ti .1 \ .1 ri.. liil h V. All vcrt.i>. 1 f>f .notion froo'. a f.-w 

.'."liin.-'.cU-rs up to tm'.s of kilonictors a.re included in these profiles. 

To dc:i'icn spat ^ \'cliicles v/ilii sulfi<dvoit slructui'a! Uaic. control systen 
; >.i ! ■ 1 1 iL y L.) Wilasliintl tin." sl(>cha;‘.lLc 'horci fiUc ti^ the f|i;l.tiiofi 

-.^l.n:cturc of these profiles (wind gusts;, the design engineer requires 
L model of wind gvxsts which captures the stochastic character of tne 
wind oroiilc and yet, at the same time, permits ease of application. 
Ti;e scales of motion associated with vertical wavelength A < 2000 m 
are responsible for the high-frequency excitations of the structural 
a'-.c; eontroi system modes of response of space vehicles [l]. The 
gusts associated with these small-scale wind profile features are the 
subject of this paper; however, before proceeding further, let us 
consider what is available to the engineer in order to focus on the 
contribution presented herein. 

Various types of gust "inputs" are available to the aerospace 
cngirxcoring comnounity for vehicle design. Included in these inputs 
are discrete gusts, spectra, and samples of v/ind profiles [2]. 

The discrete gust models generally consist of a gust shape 
(ramu builduo and decay, 1 -cosine, etc. ) with a gust amplitude 
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. I- .ill iv- . .'u.'lV'.t.-U r.ts of (!ic I’e^uUinj: tliffe ronllal equcilio.is 

1 ."o a ;; s i -I :• V' c! valuoy appropi'iatc iov a specific point in tiriic on -ho 
,'i. 1 \''c»iiclo ii'ajcctony, nncl a rosponst: unnl/sis is conducted 
iocnlly tibout iliis sclccit;ci tinic point [•!]• Pbe resulting linear, 
constU'.iT; cccfficienl, differential equations arc i ourierr transiormcc, 
resulting in a system response function from wliich, in conjunction 
wlili I'iie design gust spectrum, tnc vehicle rc^spOjise speci.ra arc 
obtained. These results can then be used to make statistical design 
judgments. To do this, a Gaussian gust process is assumed. To 
anoly the design spectrum to nonlinear vehicles with nonconstant 
coefficients, the spectrum can be factored and back-iourier transforr.»ed 
to ibe vertical coordinate domain to yield a liitcr equation which ho.s 
C‘ ui.ssian wliito noise as a lui’cing ivinction. 1 ho illli-i prov.ee s a 
series of random functions that liavc second moment statistics which 
resemble the design gust spectrum, and these random functions arc 
used as a simulation of the design gust environment. These simula.ee 
gust snacc histories are then used to "drive" the nonlinear v'ehiclc 
ecuations. A major difficulty here is that the resulting gust profiles 
do not exhibit certain vertical nonhomogeneous features. More will be 
said about this point later. 
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iiii" \ '.uli, abi.' 1’ U' i^U;^0, 8.ii;iCU:riL Iti j,t;. ^ri;'‘ 


., ^ li ;! I . V.' ' - .1 li wi'.iii ]> !’>. j ; i !(.' r^fxijjiv. . 

A critic ibin ul‘ liu; lipcctrai ir.cdcls discussed prt vLt>u.../ 

iiv . b. sumption of sluti:. i iual ■.u>:'i-iot'cne Ity alon},' the vcrticai. It 
i- V. oU -V siubl 'ccl fad Lliul ilu- v.’bifl siructnre alonj;; the vertical 
..i ‘.ut; t ropo spboi'C and lower stratospnere is not statistically hoinoc^ooeou.' . 
I'otii oxperimer.tal and thc.o rcticai evidence indicates this nonhoiViO- 
penoous character. For exaniple, on the theoretical side of the coin 
ihe basic gravity-wave theory of Hines [8] shows that the Fourier 
amplitudv?s of gravity-wa ve perturbations increase with altitude z as 
e:-;p (z/2K), where II is the atmospheric density scale height for an isothermal 
a,tmoophere. On the experimental side, extensive wind profile data 
r.re a.v?.ilable which show nonhomogeneous behavior. For example, 
the eddy structure of the wind profile in the straitosphere and troposphere 
is marlcedly different in the sense that the amplitudes of the gusts are 
;;cncrally larger above the tropopause [9]. Honhomogenoity of the wind 
profile gust statistics, and hence spectra, should be the normal state 
of affairs in view of the obvious vertical nonhomoueneous structure <}f 
c-tmes phe ric thermodynamic variables and associated synoptic (weather 
map) scale wind fields. Accordingly, the goal of the research reported 
herein is to provide a nonhomogeneous stochastic model of the 
horizontal gusts on the vertical profile of the wind which have wavelengths 


A < 2000 m. 
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1 


i.s to mode: the correlation function. This will, by necessity, ir.v.jh, e 
a iro.nsformation of the vertically ncuihomo^cr.eous ranfioin zonal and 
r.ieridional gust processes u(z) aiid v(z) Into the homogeneous (\vc hopc.i 
"recesses u(z)/a (z) = f=[z/L(z)] and v(z)/o,. (z) = ;[z/L{z)], where a (/.) 

ti ^ 

and c i'z) arc the zonal and meridional gust standard deviations at 

V 

height z. The quantity L(z) is a length .scale which is defined such that 
upon ;-;.retching transformation of the gust altitude coordinate, the 
ra.ndoni processes t(t) and ^(t) are statistically homogeneous in the 
nondimens'. -lal coordinate t = z/Li(z). A spectrum function shall be 
derived from the correlation functions of £(t) and ^(t) in terms of a 
iicndiniensional wave num.ber. Factorization of this spectrum function 
will then lead to a set of first-order state differential equations for 
^(t) and b(t) driven by white noise random processes. These state 
ecuations facilitate the sunulation of the homogeneous processes £(t) 
and .,(t) which, in turn, can be transformed to simulations of the 
ncnhoniogeneous random processes u(z) and v(z). The sirr-iulations 
of u(z; and v(z) can then be used for space vehicle design studies. 
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'14 \ c* ^ 


To Ccvoior* a .^lut M* r:UHl.r! lioM;'.onlJ.l v/iiici ;^.uhts .suitaLl 


■ ..• i..-..ilaiioii oi* V. ind proxies along ihc vertical for gust v/avclur.-t ns 
A < 20v t; iw a sample 01 cletaileu wind profile obscrvalions measured 
... intervals with ih. I'PS-l 6/ Jlmsphcrc v.'ind measurement 

-y. torn [6] was selected lor analysis. Pacii wind profile observation 
c : :'.s Lsted of a zonal and a mcridioiial wind proiile in the 1 to 18 nii. 
aliiivtde hand. The total sarriplo consisted of approximately 150 wind 
orD-i'^; observations each for the months of January, February, and 
: .rch. The zonal and meridional wind profiles for each observation 

v\ - re digitally filtered with Martin-Gr«.iim filters specially developed 
by DeMandol and Krivo [lO] for detailed wind profile analysis with a 
cutoff at X -- iiOOO m cliaracterized by a very sharp "roll-off, " so that 
^ .-.sentially all information associated with X - 2000 m was retained 
Nvith a minimum of extrancovis information associated with X > 2000 m 
retained in the gust profiles. The "DC components" ol the gust p_-* es 
wort removed by averaging the ensemble of zonal and meridionaL. 
liltercd profiles separately at each altitude and subtracting the resulting 
ensemble average zonal and meridional profiles irom each of the 
respective profiles in the ensemble. The DC components reveaicd no 
detectable nonzero trend along the vertical and were statistically equal 
to zero. Wc .denote the short wavelength (X 2000 m) zonal and 
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7'->- cmoirica: two -no int covuriuuce functions 
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(i) 


cleterinincd for z iiicrcinented at lOOO-m intervals and z incre- 
r 

’iicntcd at 25-m intervals rel.ativc to altitude z .or z < z . The' 

r r 

angular brackets in Eqs. (1) and (2) denote ensemble averages, and 

ci^;{z) = <u"(z)> , a^“ = <v“{z)\ . (3) 

Tigeros 1 and 2 contain the exiper iaaental estimates of and versus 

.’.nndimensional space lag (z^ - z)/L.(z^_). The length scale L(z^J is the 

value of z - z associated with the first zero of R. The experimental 
r 

values of and L are given in Figures 3 and 4. The data in Figures 
I tlirough 4 appear to show that the statistics of u and v arc equal, so 
that in the analysis that follows we shall pool these statistics. This is 
physically reasonable. The following formulae summarize the a and L. 
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. (•/,> s 1. ’>077 ni lii-v’ * . ' VI j 

‘:r 0. 3-i6 o.xpa.'lS X lO"*/.) . *"• 

L(;;> “ 310 + 0.0139-s, < VI ^>0 ^5 

*i2B /. u 0160 m 

^.-hcve all units arc in the MKS system. 

AUTOCOVARL^NCE AN3 SFSCTRAL FUNCTIONS 

It appears that the scaling of u(a) and v(a) indicated in Eqs. (1) 


aiid ^3) and the nondime 


tnsionaliaation of (r. - n) by Jiviaion with 


U,. ) result in a rcasonahh- collapse of the ovarinnee data and. thu-.. 
p.oidde a basis for the developntcnt of a statistical "UW. or model 

4.1 4 *t US cxnrcss the nondimcnsional 
o: detailed wind profile gusli., het us expre 

lag as 
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■ '! ^/j, l.hf lCLfi;.U: i-’i wr;;iTi ^ i 

i; M.’icI ‘ r ;t.>i ,(1 [il]. Jr. our pjirlicu’ar 
x.(l) arc linearly related to f,(t) and C(t) 


0 5 


; 


« 


where v(t) Ls either |(t) or ^ft). Auer a nu:nbcr of jnathcmatical atepci, 

it i.-; easily \erified that the coefficients a.., cl., and e. are given by 

11 1 1 



The state space system given by Eq. (14) will yield the system transfer 
function, Eq. (15). T.ne flowgraph for the state space representation of 
the gusts is given in Figure 6. This system is stable and will serve 
as a basis for developing digital filters for the simulation of zonal and 
meridional gust i^rofucs. 
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to within first order in T 


'. 2 - 


r 


Thus, the quantities ih (T) and A (T) are known functions of T 
and the system, parameters D and B, and Eq. (20) thus permits the 
(‘ijiitaal simulati(.>n of x.(K); to obtain the digital gust processes wc 
merely use Eq. (15), so that 

y(K) = e.x.(K) . • _ (25i 


DISCUSSION OF RESULTS 

To apply the above model one first generates the discrete Gaussia.n 
noise process w(K). Substitution of the noise process into the recursive 
filter given by Eq. (20) and application of Eq. (24) yields the nondimens iontil 
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rcr. V,..;- d it i ^'C-d rit ii'it-'rv-.J.};: >lon.^ tli<’ verticil, -tiVi t.'ic 

jjvo'.'O .-i :<inc procedure ubt:d to e’erive the wind proiiles irom t:-c 
‘ I'S-lo trickir.;^ d;Ha were do.si|>,ned ao that nearly all ( ~00 perccnii 
oi 5 'Oi. ctval energy of the gust process existed at wavclengtn 

- 100 vn with zero spectral energy at X = 50 m Inis means -he 

l,vv:;;ist wav'clencth is 50 m, which corresponds to a vertical oigiti- 
..lloa ii.trr\ il .Xr, '' 1’.5 :n !n /.-rpa. e. Since llie typical valiu- o: 

L, = 400 m (see Figure 5), then the digitational interval of the input 
nci; c process in t-space should be no larger than T — Az/L. — 0. Oo. 
Fi'ture 5 gives the results of a 1000-point simulation with T = 0. 06 
end is iit guod agrceincnt with the desired autocorrelation. 


CONCLUDING COMMENTS 

The statistical model reported in this paper was developed by 
assuming the gust process could be rendered nondimens ional by scaling 
it v.'ith the gust standard deviation a(z) and nondimensionalizing the 
altitude with a length scale L(z). However, classical theory of con- 
tinuous media and dimensional analysis concepts demands that 
additional parameters exist by which to scale a(z) and L(z). This is 
a crucial point because to apply the model to locations other tnan the 
Kennedy Space Center, the param.eters a(z) and L(z) must be determinec 
Nondimensionai functions for cr(z) and L(z) would permit this objective 
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